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We study the superconducting instabilities of singlet and triplet pairing in a two-dimensional 
Hubbard model on the basis of the third-order perturbation theory (TOPT). We investigate 
the effect of the vertex correction that is given by TOPT, comparing with the study with only 
the second-order effective interaction. In our results, a stable p-wave pairing state spreads from 
low to intermediate electron density. A d^2_y2-waYe pairing is dominant for the high density 
near a half-filling. It is shown that the vertex correction plays an essential role in making the 
p-wave pairing state dominant. 
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§1. Introduction 

The superconductivity originating from the electron correlation has been investigated in a quasi- 
two-dimensional system. As a simple model describing the system, a Hubbard model is used in 
many studies^^^^ for a pairing instability. The instability is often studied with a second-order 
effective interaction with respect to the Coulomb interaction. By these previous studies, it is well 
known that the dominant state is a (i^2_y2-wave pairing state near a half-filling (the density per 
spin n=0.5). The result is similar to that of the cuprate superconductors. 

The pairing instability has been analyzed also for the electron density lower than the half-filling. 
Several different consequences are suggested in some analyses. Kondo& and Chubukov et al.^'^^ 
studied about the pairing instability with the second-order effective interaction. They investigated 
the Hubbard model based on the weak coupling theory at the zero temperature. They indicated the 
advantage of a d^y- and (i^2_j^2-wave pairing symmetry. Hlubina& reported that a p-wave pairing 
state is possible in the low density at finite temperatures. He indicated that the p-, dxy- and g-wave 
pairing states are located in the lower density. 

On the other hand, Chubuko^ii considered an influence of the third-order interaction on the 
p-wave pairing at the zero temperature. They studied the problem on the basis of the Fermi 
gas model. Their results of the two-dimensional system are different from those obtained by the 
second-order interaction, especially for the p-wave pairing. They concluded that the p-wave pairing 
is induced only by the third-order interaction, and the p-wave pairing is not encouraged by the 
second-order interaction possessing a very weak dependence on the wave number. 

Nomura and Yamada0^ made also an analysis for the p-wave pairing with the third-order per- 
turbation theory (TOPT).^^~^^) They investigated a two-dimensional Hubbard model at finite 
temperatures for the density higher than the half-filling. Their conclusion is similar to the result 
which ChubukovB obtained about the Fermi gas. They concluded that the important factor for 
p-wave pairing is the wave number dependence of the third-order effective interaction, where the 
third-order interaction for triplet pairing arises from the vertex correction. 

From these results, we consider that the third-order interaction plays an important role for the 
p-wave pairing far from the half-filling in contrast with the analysis based on the second-order 
interaction. We study the pairing instability at finite temperatures for the density lower than the 
half-filling. In this paper, we investigate the problem for the two-dimensional Hubbard model on the 
basis of TOPT. We obtain a phase diagram where the parameters are the hopping integral and the 
electron density. The phase diagram is compared with Hlubina'si^ and Chubukov et aZ.'s^'^^ ones. 
We analyze the role of the third-order interaction together with that of the second-order interaction. 
Moreover, we actually solve the Eliashberg equation to obtain the transition temperature Tc of the 
superconductivity. 

Before starting our study, we should discuss the magnetic instability. In quasi-two-dimensional 
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systems, the magnetic instability is considered to be much reduced. Actuahy there exist supercon- 
ducting ground states in the quasi-two-dimensional systems such as cuprates, K-type BEDT-TTF 
and Sr2Ru04. Moreover we can expect that the superconducting fluctuations also suppress the mag- 
netic instability. However, it is difficult to determine the actual ground state by using theoretical 
calculations, since we cannot avoid using different approximations for magnetic and superconduct- 
ing states, respectively. In this paper, we confine ourselves to the superconducting states in the 
discussion on the dominant phases. 

§2. Formulation 

The two-dimensional repulsive Hubbard Hamiltonian for a square lattice is given by 

'H = -ti Yl 4,aCj> + *2 Yl 4,aCk,a + Uj2ni,^ni,i, (2.1) 

<i,j>,(T <i,k>,a i 

where a is the spin index, < i,j > indicates taking summation over the nearest-neighbor sites 
and < i,k > over the next-ncarest-neighbor sites. We obtain the energy dispersion from the 
non-interacting part in eq. (2.1); 

Ek = —2ti{coskx + cosky) + At2COs{kx)cos{ky) , (2.2) 

where we take ti=1.0 and —0.5 < t2 < 0.5. By using the above dispersion, wc obtain the bare 
Green's function given by Go{k,en) = ' where = 7rT(2n + 1) is the Matsubara fre- 

quency and fi is the chemical potential. The particle number (density) per spin n is given by 

n = ■^^k,nGo{k, e„). 

An effective pairing interaction in this study is given by the perturbation expansion up to the 
third-order term with respect to the on-site Coulomb interaction. Wc calculate the effective inter- 
actions for the singlet and the triplet states, respectively. The effective interaction is divided into 
two parts, 

VTOPT(g, k) = V-RPAiq, k) + Vvertex(q', k). (2.3) 

The RPA-like term Vrpa includes the term given by the Random Phase Approximation (RPA) and 
Vvertex is the vertex correction. The RPA-like term reflects the nature of spin fluctuations. The 
vertex correction term originates from the electron correlation U other than the spin fluctuations. 
For the singlet pairing, the RPA-part and the vertex correction part are given, respectively, by 

V"S'*(9, k)=U + U\^{q -k)+ 2U\l{q - k), (2.4) 

^Vertex* (9, k) = 2^U'[Y: Go{q -k + k')x {xoiq + k') - Mq + k'))Goik% (2.5) 

k' 

where k indicates k = {k, tOn)- The bare susceptibility Xoiq) 4'oiq) ^^^^ deflned respectively by 

Xo{q) = E ^^0(^)^0(9 + k), (2.6) 

k 
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Mq) = E Go(k)Go{q - k). (2.7) 

k 

For the singlet pairing, the Coulomb interaction U connects only the electrons which have the 
opposite spins. The diagrams for the pairing interaction is shown in Fig. 1. The two external lines 
have the opposite spins. 

For the triplet pairing, the RPA-like term and the vertex correction correspond to the second- 
order and the third-order terms, respectively. 

V^Jt\q,k) = -U\o{q-k), (2.8) 

VyJ^{q,k)=2^U'[Y,Go{q-k + k') X (^.{q + k') + Mq + k'))Go{k')]. (2.9) 

k' 

As an effective interaction for the triplet pairing, the diagrams are shown in Fig. 2. The two 
external lines have the parallel spins. 

An anomalous self-energy is written by the effective interaction V{q, k) and an anomalous Green 
function F{k), as S^(g) = —j^Ti^F {k)V {q, k). At the transition temperature T = T^, the value of 
the anomalous self-energy is small and we linearize the Eliashberg equation with respect to F 
and S^, as F(fc)^ = |Go(/c)|^SA(fe)^- From these formulae, we obtain the following equation for 
the anomalous self-energy; 

^4(9) = ~Y.^{qMGom^^\{k). (2.10) 
k 

This equation is the linearized Eliashberg equation which is an eigenvalue equation with an eigen- 
value A and an eigenvector S^. V{q,k) is given by (2.3), (2.4), (2.5), (2.8) and (2.9). We solve the 
linearized Eliashberg equation on the assumption that S]^ has the pairing symmetry represented by 



p-wave; sin(ky), 

c?2,2_j^2-wave; cos(kx) — cos(ky), 

dj^y-wave; sin(kx)sin(ky), 

g'- wave ; sin (kx ) sin (ky ) (cos (kx ) — cos (ky ) ) . 



When the eigenvalue calculated from eq. (2.10) reaches unity, the superconducting state is 
realized. The most dominant pairing symmetry has the largest value of the eigenvalues among 
different symmetries. We solve the equation to obtain the dominant state and determine at 
which the eigenvalue A equals unity. 

Here, we add a comment on neglecting the normal self-energy Sji(g, oj) in our calculation. When 
we calculate Sn(g, w) with TOPT, Sji(g, a;) changes rather drastically for the density far from the 
half-filling. The change is due to the large asymmetry with respect to electron-hole properties. In 
the low density case where the numbers of electron and hole are very different, the real part of 
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the normal self-energy IieTi^{q,uj) does not show the smooth Fermi liquid behavior. The wrong 
behavior of Sn(g*, is due to the fact that Sn(g*, is confined to a finite order. The behavior 
oscillates with respect to the order of U, owing to a weak convergence. 

To avoid the difficulty due to the weak convergence, we use the following method treating the mass 
renormalization z. The method is discussed in the theory of heavy fermions.© By this method, 
we treat separately the frequency and the momentum dependencies of ReSn(g, w). The frequency 
dependence of ReSn(g, lij) is included by z in the starting Hamiltonian, where z is determined so 
as to give a correct effective mass enhanced by z~^. Then, the momentum dependence is included 
in the perturbation calculation. 

The detail of the method is the following. The bandwidth W of quasi-particles is renormalized by 
z. The interaction between the quasi-particles F is renormalized as zTz = zU, where F is enhanced 
as F = [/ /z. As a result, W and F are renormalized by multiplying z. Therefore, the Hamiltonian 
(2.1) can be rewritten with ti, t2 and U including the renormalization factor z. (We don't exhibit 
the rewritten Hamiltonian explicitly here.) Thus, we can include the frequency dependence of 
ReSn(^, uj) in the starting Hamiltonian. In the procedure for the frequency dependence, we neglect 
the momentum dependence of ReSn(^, w). On the other hand, the momentum dependence of 
the effective interaction is essential in realizing the anisotropic superconductivity. The momentum 
dependence is included by solving the Eliashberg equation obtained by the perturbation calculation 
with respect to U, as shown in this paper. The separation between frequency and momentum 
dependencies does not change the essential physics and the procedure of mass renormalization 
becomes easy. 

§3. Numerical Calculation and Results 

We divide the first-Brillouin-zone into 256x256 momentum meshes and take A^f = 1024 for 
Matsubara frequency. To make a reliable calculation, the region of Matubara frequency a;„ should 
cover the bandwidth W=8. The region is covered with the condition; |u;n| < W. To meet the 
condition, our calculation is confined to the temperature region T > 0.002. In addition to the 
condition, the value of U is confined to U < W. 

3.1 Superconductiong phase diagram and xoiQ,^n = 0) 

In Fig. 3, we show the phase diagram for the dominant pairing state in the plane of the second- 
nearest-neighbor hopping integral t2 and the electron density n. The regions of parameters are 
given by 0.075 < n < 0.4 and —0.5 < t2 < 0.5. The temperature T equals 0.008 and the Coulomb 
interaction is U =6.0. In Fig. 4, the dependence of the bare susceptibility XoiQ,^n = 0) on n is 
shown for a quarter first-Brillouin-zone. 

For the high density (n > 0.3) near half-filling in the phase diagram, the (i^2_y2-wave pairing 
state is dominant due to antiferromagnetic fluctuations which have a peak of susceptibility xo at 
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(tt, tt) in the momentum space. The situation is the same as that of cuprates corresponding to a 
nearly half-filled system. The ^?-wave pairing state is realized for the density from n = 0.075 to 
0.303. The stable p-wave pairing state spreads from low to intermediate density (0.075 < n < 0.3). 

3.2 Dependence of X and Tc on n 

We show the dependence of the eigenvalue A on n in Fig. 5. We solve eq. (2.10) with r=0.008, 
t2=0.25 and ^7=6.4. The region of n is given by 0.03 < n < 0.375. The most dominant pairing 
state has the largest value of A among various symmetries. 

Fig. 5-(a) shows the eigenvalue for eq. (2.10) with Vtopt for all terms which are given by 
TOPT. Thep-wave pairing state is dominant over the other states from low to intermediate density 
(0.075 < n < 0.3). For the high density (n > 0.3) near the half-filling, the d^2_j^2-wave pairing 
state is realized. The eigenvalues for the other singlet states (the dxy- and the g-wave pairing) 
never reach unity which gives the transition temperature Tc. 

Fig. 5-(b) shows that A obtained from eq. (2.10) with Vrpa given by the RPA-like term, which 
includes the RPA-like term reflecting the spin fluctuations. The (i^2_y2-wave pairing state is realized 
for the high density and the dxy-'wave pairing state is advantageous from intermediate to low density, 
respectively. In the case of the p-wme pairing state, A given by RPA-like term is smaller than one 
given by TOPT. (The RPA-like term for the triplet includes only the second-order term.) For the 
singlet pairings, A obtained from the RPA-like term are larger than one of TOPT. From a unity 
hue of A in Fig. 5-(a) for TOPT, we find Tc=0.008 for the p-wave pairing state at ri=0.075 and for 
the dx2_y2-ws,ve pairing at n=0.289. 

In Fig. 6, we discuss the detail of the eigenvalue calculation and make clear each role of the 
RPA-like term Vrpa, the vertex correction VVertex and TOPT Vtopt- The p-wave pairing state has 
the largest eigenvahic for the case possessing only the vertex correction in Fig. 6- (a). (The vertex 
correction for the triplet includes only the third-order term.) The eigenvalue A for the RPA-like 
term is the smallest. The eigenvalue A given by TOPT became smaller than one given by the 
vertex correction. The large eigenvalue obtained from the vertex term is suppressed by the RPA- 
like term. This is because that a scattering amplitude of ^^"tex* opposite phase to that of 
y^ipict^ About the singlet states in Fig. 6-(b), (c) and (d), the RPA-like term promotes the d- and 
g-wave pairing. Thus, these singlet states are encouraged by the spin fluctuations. For example, 
the antiferromagnetic fluctuations are the main origin for an attractive force of the d^2_y2-wa.ve 
pairing. In contrast to the RPA-like term, the vertex correction suppresses the eigenvalue and 
lowers the transition temperature. 

3.3 Transition temperature Tc 

We show the superconducting transition temperature Tc in Fig. 7. In Fig. 7- (a), we show Tc 
which is calculated with TOPT. The unit of energy is the hopping transfer ti. Tc for the p-wave 
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pairing state is lower than Tf. = 0.07. Tf. for the d^2_y2-wave pairing state is higher than T^. = 0.03 
for the high density (n > 0.3) near the half-fiUing. In Fig. 7-(b), we compare Tc which are obtained 
by TOPT, the RPA-Hke term and the vertex term, respectively. For the p-wave pairing symmetry, 
Tc obtained by only the vertex term is higher than that of TOPT by 0.03. On the other hand, Tc 
given by the RPA-like term for the d^2_y2-waYe pairing symmetry is higher than that of TOPT by 
0.2. These fact means that the p-wave pairing is induced by the vertex correction and dj.2_j,2-wave 
pairing is induced by the RPA-term. 

3.4 Dependence of X on U 

In Fig. 8, 9 and 10, we show A at U= 2.0 and 4.0 to study a dependence of A on U. The 
hopping integrals are fixed as ^2=0. 25 and 0.4. Fig.8 exhibits A obtained by TOPT for the p- and 
dx2^y2-wave pairing. About the result for t2=0.25, the dominant state at U = 6.4 is the p-wave 
pairing for about n < 0.3 and the dx2_y2-waye pairing in the high density (n > 0.3) in Fig. 5-(a). 
The situation docs not vary when the value of U changes to 2.0 and 4.0 in Fig. 8-(a). For t2=0.4 
(Fig. 8-(b)), the situation docs not also vary for the value of U. Thus, the value of U does not 
change the dominant state for the various density in the calculation for TOPT. 

Fig.9 shows A obtained by RPA-hke term. At U=2.0 and 4.0 for both ^2=0. 25 and 0.4, the 
advantageous pairing state is the dxy-pairing state from low to intermediate density (0.075 < n < 
0.3) and the (i3,2_j^2 -pairing state in the high density. However, the p-wave pairing is dominant for 
the very low density (ra < 0.075) in the case of i2=0.25 for U=2.0 and 4.0. For [7=6.4 in Fig. 5-(b), 
the (ia;y-pairing is dominant for n < 0.275. The p-wave pairing is not advantageous for the very low 
density. When U decreases from 4.0 to 2.0, the density for the p-wave pairing extends from 0.03 
to 0.075. In the very low density for small U, we think that the p-wave pairing becomes dominant 
owing to the effect of the ferromagnetic spin fluctuation (paramagnon) or too small eigenvalue for 
the comparison. 

Fig. 10 exhibits the [/-dependence of the eigenvalue A obtained by including the all (TOPT), 
vertex and RPA-like terms for the p-wave pairing state. Similarly to Fig. 6-(a) for [7=6.4, the 
results in Fig. 10 show that the vertex correction mainly promotes the p-wave pairing state. The 
RPA-like effective interaction suppresses A for the p-wave pairing. This situation does not change, 
when [7 equals 2.0 and 4.0 in the cases of both t2=0.25 and 0.4. 

The results of Figs .8, 9 and 10 indicate that the situation of the dominant pairing state does 
not change for the value of [7, except for the very low density in the case where U is smaller than 
the half of the bandwidth. 

3. 5 Wave number dependence of V 

The wave number dependence of the effective interaction V in the momentum space is shown in 
Fig. 11. The wave number dependence brings an electron scattering near the Fermi surface. The 
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scattering causes the attractive force between electrons on the Fermi surface to reahze the pairing 
state. The point k = kp on the Fermi surface and q denote the initial state and final state of the 
scattering, respectively. The lighter color shows a region of the stronger effective interaction which 
gives the frequent scattering from kp- 

When the (i^2_y2-wave pairing state is advantageous, the singlet effective interaction is shown 
in Fig. 11- (a). The RPA-like interaction V^p^^''* reflects the antiferromagnetic fluctuations and 
gives the scattering leading to the dx2_y2-wave pairing state near the Fermi surface. The vertex 
correction V^grtex* makes the scattering weak for the ^^.'^-yZ-wave pairing. When the p-wave pairing 
state is dominant, the triplet effective interaction is shown in Fig. ll-(b). The vertex correction 
^^rtex* iiiduces the scattering leading to the p-wave pairing near the Fermi surface. The RPA-like 
interaction make the scattering weak. Thus, the RPA-like interaction never gives the strong 

scattering enough to lead the triplet pairing. 

§4. Summary and Conclusions 

We have used the third-order perturbation theory for the two-dimensional Hubbard model and 
studied the dominant pairing state for the superconductivity. We compare our result with the 
previous one that includes only the second-order effective interaction. 

In contrast with the case including only the second-order interaction, the region of the p-wave 
pairing state is extended broadly from low to intermediate density(0.075 < n < 0.303). Thus, the 
introduction of the third-order terms gives the important influence on the p-wave pairing. The 
dxy- and g-wave pairing has been suggested in Hlubina's papeii^ based on only the second-order 
calculation. However, these singlet states do not appear as the advantageous state in our result by 
TOPT. 

We discuss more detailed consequences in the followings. In the triplet state, the vertex correction 
brings mainly the advantage for the p-wave pairing. The vertex correction has the wave number 
dependence. The wave number dependence causes the scattering near the Fermi surface, which 
induces the attractive force for the p-wave pairing. On the other hand, the RPA-like term suppresses 
the p-wave pairing. The spin fluctuations lower the eigenvalue of the p-wave pairing. Thus, the 
ferromagnetic spin fluctuations do not play a role in realizing the p-wave pairing. The situation 
from low to intermediate density is different from that in the very low density for small U. In the 
very low density, the ferromagnetic spin fluctuations (parramagnon) encourage the p-wave pairing. 
However, the eigenvalue for the paramagnon mechanism is less than that for the vertex correction. 

In the case of the singlet state, the vertex correction discourages the singlet pairing. Thus, the 
vertex correction lowers the transition temperature. The singlet state becomes dominant by the 
RPA-like terms reflecting the spin fluctuations, which has been shown by many previous studies. 

Our result might interpret the origin of a quasi-two-dimensional superfluidity of 3jjg8,9,i6)_ 
dicate that the p-wave pairing becomes dominant owing to the paramagnon in the very low density 
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case. The results agree with the consequence about a two-dimensional superfluidity of He obtained 
by previous studies.^ However, the density corresponding to the realistic two-dimensional ^He- 
superfluidity is not the very low density. For example, the density of a realistic three-dimensional 
^He-superfluidity is the quarter-filling near the intermediate density. Therefore, we think that 
the density of the two-dimensional ^He-superfluidity is the intermediate density. Thus, the ver- 
tex mechanism might give the important effect on the superfluidity rather than that due to the 
paramagnon. Our results might suggest that the vertex correction plays the important role for the 
superfluidity of '^He as well as the ruthenate superconductor. In order to make clear the origin 
of the ^He-superfluidity, it is necessary to consider not only the paramagnon but also the wave 
number dependence of the vertex correction. 

In addition to this, we mention the effect of the lattice for the intermediate density. The influence 
of the square lattice is weak in the very low density, because the Fermi surface in the very low density 
is similar to that of the Fermi gas. On the other hand, the effect of the lattice can not be neglected 
in the case where the density is not very low. Therefore, it is necessary to compare the results 
obtained for various lattices on the basis of TOPT. 

Our result might suggest the reconsideration concerning the origin of not only two-dimensional 
^He superfluidity but also the three-dimensional one0 A part of this paper is to be published in 
proceedings of ISSP meeting® 
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Fig. 1. Diagrams for the effective interaction of singlet pairing within the third-order perturbation with respect to 
U. The solid line is the bare Green's function Go. The broken line is the Coulomb interaction U. The broken line of 
U connects only solid lines possessing opposite spins. The two external lines have the opposite spins. The effective 
interaction is divided into the RPA-like part and the vertex correction, the latter beginning from the third-order 
terms. 
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Fig. 2. Diagrams for the effective interaction of the triplet pairing within the third-order perturbation. The two 

external lines have the parallel spins. The RPA-like part and the vertex correction of the pairing interaction are 
given by only the second-order and the third-order terms, respectively. 
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Fig. 3. A phase diagram for the dominant pairing symmetry. The diagram has the parameters of the second-nearest- 
neighbor hopping integral t2 and the density n. t2 is from -0.5 to 0.5 and n is from 0.075 to 0.4. (A half-filUng 
corresponds to the density n=0.5.) The temperature T equals 0.008 and the Coulomb interaction U equals 6.0. 
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Fig. 4. A bare susceptibility xo{<l,^n = 0). The dependence of the susceptibility xo(5, = 0) on the density n in 
a quarter first- Brillouin-zone. The parameter sets are similar to those of Fig. 1. 
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Fig. 5. The dependence of the eigenvalue A on the density n. (a) is calculated with all the terms given by the third- 
order perturbation (TOPT). (b) is obtained from only the RPA-likc term. The parameters are the temperature 
T=0.008, the hopping integral t2=0.25, and the Coulomb interaction U=6A. The region of the density n is from 
0.03 to 0.375. 
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Fig. 6. The dependence of the eigenvahic A on the density n for (a) p-, (b) d^2_y2-, (c) dxy- and (d) g-wave pairing 
state. The parameter sets are similar to those of Fig. 5. 
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Fig. 7. The superconducting transition temperature Tc. (a) Tc which is calculated with all terms (TOPT) given 
by the third-order perturbation theory, (b) The comparison of Tc obtained with all terms (TOPT), the RPA-like 
term (RPA) and the vertex term (Vertex). The unit of energy is the hopping transfer integral ti. 
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Fig. 8. The {/-dependence of the eigenvalue A obtained 
theory for the p-, d^2_y2- and dxy-vfave pairing states, 
of U are 2.0 and 4.0. The temperature T equals 0.008. 



from all terms (TOPT) within the third-order perturbation 
The hopping integral is fixed as t2=0.25 and 0.4. A values 
The region of the density n is from 0.03 to 0.375. 
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Fig. 9. The {/-dependence of the eigenvalue A obtained from the RPA-hke term for the p-, d^2_y2- and da^^-wave 
pairing states. The hopping integrals are fixed as t2=0.25 and 0.4. The values of U are 2.0 and 4.0. T equals 0.008. 
The region of the density n is from 0.03 to 0.375. 
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Fig. 10. The [/-dependence of the eigenvalue A obtained from the all (TOPT), vortex and RPA-like terms for the 
j3-wave pairing state. The hopping integral is fixed as ^2=0. 25 and 0.4. The values of U are 2.0 and 4.0. T equals 
0.008. The region of the density n is from 0.03 to 0.375. 
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Fig. 11. The dependence of the singlet- and triplet-effective interaction V on the wave number q for all the terms 
(TOPT) Vtopt, the RPA-like term Vrpa and the vertex term Wertex on the momentum space q. The point k = kp 
fixed on the Fermi surface and q are an initial and a final state of the scattering. A lighter color shows a region of the 
stronger effective interaction, (a) The singlet effective interaction in the d3.2_j,2-wave pairing state. The parameter 
sets are T=0.008, U=6A, ^2=0. 25 and n=0.225. (b) The triplet effective interaction in the p-wave pairing state. 
The parameter sets are T=0.008, U=6A, t2=0.25 and n=0.375. 
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